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ABSTRACT 

We present the results of high angular resolution adaptive optics (AO) near- 
infrared {JHKg) observations of the deeply embedded massive cluster RCW 38 
using NACO on the VLT. Narrowband AO observations centered at wavelengths 
of 1.28 /im, 2.12 /im, and 2.17 fim were also obtained. The area covered by 
these observations is about 0.5 pc square, centered on the O star RCW 38 IRS2. 
We use the JHKg colors to identify young stars with infrared excess in this 
region. Through a detailed comparison to a nearby control field, we find that 
most of the 337 stars detected in all three infrared bands are cluster members 
(~317), with essentially no contamination due to background (hkely due to the 
high cluster extinction of Ay ~ 15) or foreground sources. Five sources with 3 
band detections have colors suggestive of deeply embedded protostars, while 53 
sources are detected at Kg only; their spatial distribution with respect to the 
extinction suggests they are highly reddened cluster members but their evolu- 
tionary status is unclear. Detectable i^^-band excess is found toward 29 ± 3% 
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of the stars. For comparison to a similar area of Orion Nebula Cluster observed 
in the near-infrared, mass and extinction cuts are applied, and the excess frac- 
tions redetermined. The resulting excesses are then 25 ± 5% for RCW 38, and 
42 ± 8% for Orion. RCW 38 1RS2 is shown to be a massive star binary with a 
projected separation of ~ 500 AU. Two regions of molecular hydrogen emission 
are revealed through the 2.12 /im imaging. One of these shows a morphology 
suggestive of a protostellar jet, and is clearly associated with a star only detected 
at H and Kg, and previously identified as a highly obscured X-ray source. Three 
spatially extended cometary-like objects, suggestive of photoevaporating disks, 
are identified, but only one is clearly directly infiuenced by RCW 38 IRS2. The 
structure of the inner core of RCW 38 is also characterized and compared to 
Orion and other clusters. A King profile provides a reasonable fit to the cluster 
radial density profile and a nearest neighbor distance analysis shows essentially 
no sub-clustering. 

Subject headings: stars: formation - infrared: stars - ISM: individual (RCW 38) 



1. Introduction 

The study of embedded clusters is vital to our understanding of star formation because 
most stars have been shown to form in large clusters with hundreds of stars (Lada & Lada 
2003, Porras et al. 2003, Carpenter 2000 and Allen et al. 2007). Of particular interest are the 
largest clusters with thousands of stars because it is in these clusters where large numbers 
of low mass stars live in close proximity to luminous 0-stars. In the Lada & Lada 2003 
tabulation of clusters within 2 kpc, there are two such larger embedded clusters with over 
1000 stars: the Orion Nebula Cluster (ONC) and RCW 38 (Wolk et al 2006). The other 
clusters within 1 kpc typically contain only a few hundred stars, making direct comparisons 
with the ONC not possible. It is even likely that our Sun was formed in such a cluster. A 
possible explanation for the abundances of many radioactive isotopes in our Solar System 
is the presence of a nearby supernova during its formation which would have been a more 
probable occurrence in a rich cluster hke the ONC or RCW 38(Meyer & Clayton 2000). 

Studying embedded star clusters not only yields important information regarding the 
formation and evolution of stars, it also yields important information regarding the formation 



^ Based on observations performed at the European Southern Observatory Very Large Telescope on Cerro 
Paranal, Chile, under program ID 70.C-0400(A). 
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and evolution of planetary systems around these stars. As first studied in the Orion Nebula 
Cluster (ONC), proplyds (PROto PLanetarY DiskS) exist as large disks of dust and gas 
that provide material for future planet formation. However, when exposed to unattenuated 
ultraviolet (UV) radiation the molecular hydrogen dissociates causing the outer portions of 
the disk to photoevaporatc (HoUcnbach et al. 1994). Since many young low-mass stars in 
other regions exist in close proximity to massive OB stars that output a significant amount 
of UV radiation, it is important to determine whether this process seriously inhibits planet 
formation around them. Low-mass stars are important because they represent the majority 
of young cluster members, which means the vast majority of planets may be formed around 
them. Thus, low-mass stars clustered around OB sources are at the same time the most likely 
to form planets and the most likely to experience photo evaporation effects. Knowledge of 
the disk frequency in a cluster as a function of the distance from the massive OB stars 
can constrain the affect of photocvaporation on disk evolution (Balog et al. 2007). Such 
knowledge will help to constrain the conditions that are necessary for those disks to evolve 
into planets. 

The HII region ROW 38, at a distance of 1.7 kpc (Beck, Fischer & Smith 1991), contains 
an embedded cluster centered around a 05 star 1RS2. RCW 38 is more embedded than the 
ONC and due to the OB stars in the center it provides a similar environment for the study of 
disk photoevaporation and evolution. Smith et al. (1999) showed that the winds from IRS2 
have created a 0.1 pc bubble that is relatively free of dust, leaving stars within the bubble 
directly exposed to UV radiation. This bubble is filled with diffuse thermal and synchrotron 
X-ray emission (Wolk et al. 2002) and the intersection of the central gas bubble with the 
surrounding dust cloud appears to be a site of active star formation (Smith et al. 1999). 

Wolk et al. (2006) observed the central region of RCW 38 at X-ray and near-infrared 
wavelengths, identifying a large cluster population of 360 X-ray sources that are likely cluster 
members over the total field down to Log[L3.]~30. Comparing this to the COUP (Chandra 
Orion Ultradeep Project) measurements of the ONC discussed in Feigelson et al. (2007), 
roughly 200 sources of their sample are more luminous than Log[L2:]~30. The number of 
sources in the COUP measurement is on the same order as Wolk et al. (2006). The somewhat 
larger numbers in the latter study may be explained by the larger physical area surveyed 
relative to the former. Since the ONC and RCW 38 have similar populations of young stars, 
this shows that both the ONC and RCW 38 are unique large clusters within 2 kpc. For the 
central region of RCW 38, Wolk et al. (2006) determined that 97% of the ~200 X-ray sources 
within 200" of IRS2 are likely cluster members. In this region they determine that the disk 
fraction is 49% from near-infrared (JHKs) excess measurements. The infrared observations 
of Wolk et al. (2006) of the central 2'.5 of the RCW 38 cluster were performed with ISAAC 
at the Very Large Telescope (VLT), with resolution insufficient to resolve close sources in 
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the most crowded region around the 05 star. We present here the results of our observations 
of the central ~1' of RCW 38 with adaptive optics near-infrared imaging with the NACO 
system on the VLT, which allow us to better separate close sources and to separate the 
source emission from the nebular background emission, providing a sample that is more 
easily compared to the nearer ONC. In §2 we describe the observations, in §3 we discuss the 
photometery and completeness limits, and in §4 we present the results on infrared excess 
sources, cluster membership, the central O star, protostellar candidates, molecular hydrogen 
emission, proplyd candidates, and the disk fraction. In §5 we analyse these results, discussing 
cluster structure and the effect of the O stars on disks around low mass stars. Our summary 
is presented in §6. 



2. Observation and Data Reduction 

2.1. NAOS-CONICA 

High angular resolution near infrared imaging of RCW 38 was performed on Febru- 
ary 20, 2003 and February 23, 2003 using NAOS-CONICA (NACO) at Yepun, one of the 
8.2 m unit telescopes of the Very Large Telescope (VLT) installed at the European South- 
ern Observatory's (ESO) site on Cerro Paranal, Chile. NAOS (Nasmyth Adaptive Optics 
System) is an adaptive optics (AO) system that assists imaging at wavelengths in the range 
0.8-5.5 /im by directing light from the telescope partially into a (either optical or infrared) 
wavefront sensor. The remaining light is directed into the CONICA infrared camera that is 
equipped with a 1024 x 1024 pixel Aladdin InSb array detector (Lenzen et al. 2003; Rousset 
etal. 2003). 

The CONICA S54 setup (pixel scale ~ 0'.'054; instantaneous field-of-view ~ 56" x 56") 
was used to take 3 sets of 14 frames each in the broad-band J, H and Kg filters, centered at 
a wavelength of 1.27/im, 1.66 /xm and 2.18 //m, respectively. In addition, 3 sets of 12 frames 
each were taken in narrow-band filters centered at the following wavelengths: 1.28 /im (band- 
width = 0.014 /zm), 2.12 /im (bandwidth = 0.022 ^m) and 2.17/im (bandwidth = 0.023 /xm), 
covering the Pa (3, H2 (1-0) S(l), and 5r7 lines, respectively. Detector integration times 
(DITs) and NDITs are given in Table [TJ The integration time per frame is DIT x NDIT. 

Wavefront sensing was performed with the IR sensor to close the AO loop on the bright 
infrared source RCW 38 IRS 2 (self-referencing) at RA = 08'^59™05"5, DEC = -47''30'39'.'4 
(J2000), and SpT = 05V (Frogel & Persson 1974; Smith etal. 1999). The airmass ranged 
from 1.15 to 1.30 during the observations. For both nights, observing conditions were excep- 
tionally good and stable, with a clear sky and optical seeing mostly in the range 0"4-0'.'6. As 
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measured on the obtained data sets themselves, the image quahty (i.e., the median FWHM 
PSF of randomly selected point sources) is found to be better than 0'.'2 for all six filter 
settings. 

To allow corrections for cosmic rays and bad pixels, a jitter pattern (random offsets 
within a jitter box of width 2") was applied, resulting in a total field of view of GC.'S x 58'.'3. 
This corresponds to linear scales of 0.5pcx0.48pc at the distance of 1.7 kpc. Twilight 
fiats, lamp fiats and dark frames were taken through the usual NACO calibration plan. For 
each individual frame, all basic steps of data reduction (i.e., fiat-fielding, sky subtraction 
and bad pixel correction) were performed using standard routines within the IRAF software 
package. For an exhaustive description of the processing of crowded field IR data we refer to 
Niirnberger & Petr-Gotzens (2002). Figure [T] shows the final data set as a combined 3-color 
image. 



2.2. SOFI Control Field 

In order to estimate the degree of background contamination, a control field (J. Ascenso, 
private communication) was imaged on June 17, 2005 in the J, H, and Ks passbands using 
the SOFI (Son of ISAAC) on the NTT (New Technology Telescope) at La Silla, Chile. The 
resolution of the image is approximately 0.29 arcsec/pixel and centered at 10'^18™01';2 and 
— 57°15'31'.'l. A total of 6 frames per passband were combined together, resulting in an 
integration time of 36 s per band. The total FOV is approximately 4' x 4'. 

3. Source Identification and Photometry 

3.1. Phot Vis Identification 

Point source detection and photometry was accomplished using Gutermuth's IDL pho- 
tometry and visualization tool. Phot Vis version 1.10 (Gutermuth et al. 2008). Phot Vis 
selects point sources from embedded regions containing significant amounts of nebulosity 
with a modification to the standard DAOphot algorithm. Phot Vis is heavily dependent on 
the FITS procedures and DAOphot version 1 port supplied by the IDL Astronomy Users 
Library (Landsman 1993). 

Gaussian distributions were assumed for all sources and sources with a full-width at 
half-maximum (FWHM) greater than 2.5 pixels (0'.'135) were discarded. The images were 
also visually inspected to remove false positives due to image artifacts. Flux measurements 
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were taken using a 3 pixel (0.162") radius and background measurements were determined 
using inner and outer sky annuli of 5 pixels (0'.'27) and 8 pixels ((y.'432) respectively. This 
procedure was used to perform relative photometry on both NAOS-CONICA data and the 
SOFI control field. Absolute photometry was not performed on the NACO narrow band 
images. 

3.2. Completeness and Detection Limits 

Differential completeness limits of 90% were calculated in all images by inserting succes- 
sively fainter fake stars with gaussian profiles and FWHM equivalent to those in the image 
and taking 90% of the false sources recovered. For this test, sources with magnitude errors 
greater than 0.2 were discarded as well as sources that deviated in brightness by 3a more 
than the original (due it being overlayed over a brighter source). Final values for the com- 
pleteness limits of the NACO data are J=19.01 mag, H—1S.13 mag, and Xs=16.97 mag. 
Detection hmits {5a) in these bands are J=20.3 mag, H—19.0 mag, and Kg—lS.O mag. The 
90% completeness limit for the SOFI control field calculated in the same manner is Kg— 15.67 
mag. J and H completeness limits were not calculated because only the Kg-hand control 
field was necessary for background contamination estimates. The saturation limits for both 
fields was approximately 11th magnitude in all 3 passbands. 

Utilizing the J and H completeness limits yields a mass detection limit of 0.12 Mq at 1 

Myr for an average extinction of Ay = 10 {Ak = 1.1; Rickc and Lebofsky 1985). The mass 
detection limit was approximated using the Baraffe 1 Myr isochrone (Baraffe et al. 2002) 
and the distance modulus of 11.15 calculated with the assumed distance of 1.7 kpc. This 
detection limit is right above the hydrogen burning limit of 0.08 Mq, so there are no brown 
dwarf candidates in the NACO sample. 

3.3. NAOS-CONICA Calibration 

The NACO WCS and zero-point flux calibrations were performed using larger FOV 
images of RCW 38 obtained with ISAAC (Infrared Spectrometer And Array Camera) on 
the VLT. We used the photometry extracted from the ISAAC data by Ascenso (2008). Due 
to the use of PSF fitting and to a better estimate of the sky contribution on the nebulous 
background, the Ascenso (2008) photometry is believed to be more reliable than the aperture 
photometry of Wolk et al. (2006). The ISAAC data was registered and calibrated using 
sources from the 2MASS Point Source Catalog in the larger ISAAC field. The SOFI control 
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field was flux and position calibrated using 2MASS sources at measurements up to magnitude 
14 in JHKg, due to the fact that the SOFI field is much deeper. All sources in all images with 
a measured uncertainty greater than 0.2 magnitudes (outside of la) have been discarded. 

In addition, because the AO feedback loop was closed using a central source for the 
NACO data, there is a significant amount of off-axis elongation in many sources in the 
images. This caused a linear trend in magnitude difference between the ISAAC and NACO 
photometric measurements as a function of radial distance from IRS2. We therefore used 
position dependent zero points given by the equations: 

J : 24.426 - 0.0176 x r^s 

H : 22.828 - 0.0247 x r^s (1) 
Ks : 21.648 - 0.0274 x r^s 

where rofr is the radial offset from IRS2 in arcseconds. 

4. Results 

4.1. Color-color and Color-magnitude Diagrams 

The infrared photometry for the 483 sources detected in at least one band with a mag- 
nitude uncertainty of < 0.2 is listed in Table [31 The number of sources detected in each of 
the three broad-bands are shown in Table [21 We detected 409 sources in H and Kg and 344 
sources in all bands with uncertainties < 0.2 magnitudes. This uncertainty cutoff results in 
average photometric errors of 0.07 mag for J and 0.06 mag for H and Ks. Inspection of 
the 344 sources detected in JHKg shows that seven sources have (J — H) values less than 
1 magnitude and deviated significantly from the bulk of the sample. Of these seven, four 
have a i^^^-band excess. Upon visual inspection, it was evident that two of these sources had 
an excess of flux detected in J due to artifact contamination in that image (where residual 
images of the central source IRS2 can be seen in multiple places in the central portion of 
the image). Along the same lines, four of the other sources were very near IRS2 or partially 
obscured by a much larger, brighter star, which explains the excess flux detected in J. The 
final source, that remains unexplained by the previous explanations, is likely a foreground 
source (especially due to the fact that it does not show K^-band excess) and thus, is not 
a cluster member. Of the seven sources, six have incorrect photometry, and are tagged in 
Table [3] with an 'x'. These sources have been omitted from the final data set, leaving the 
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total number of JHKg sources at 337. It is this data set of 337 sources that is used for 
the rest of the analysis in this paper. The color-color diagram for these sources is shown in 
Figure [2] and the color-magnitude diagram in Figure Ei 

4.2. The IRS2 O Star Binary 

IRS2, the central source of RCW 38, has been previously thought to be a single source, 
but with the increase in resolution enabled by the AO system, it is possible to show that 
it is in fact a binary system. The separation of the two sources can be confirmed in all 
three broad-bands through visual inspection, but due to the brightness of the two sources, 
the detector is saturated in this region for the broadband JHKg filters. Therefore, the 
narrowband images (2.17//m and 2.12^m) have been used for a relative photometric analysis 
of IRS2. The ratio of the detected flux from the brightest pixel in each source shows the 
dimmer source to have ~ 94% of the flux of the brighter source (89% in the case of the 
1.28/im image). Comparing this ratio to the ratio of Lyman continuum photons given for 
various ZAMS spectral types in Panagia (1973; Martins et al. 2005), shows that both stars 
in the binary should have the same spectral type. Figure H] shows magnified images of IRS2 
with contours in these two bands. From the narrowband images, the projected separation 
between the two sources is 0'.'27 which is 459 AU at 1.7 kpc. 

Smith et al. (1999) address the question of the effective temperature and hence spectral 
type of IRS2 by treating it as a single source and fitting nebular spectroscopic lines to models 
of HII regions using stars with Kurucz stellar atmosphere models, from which they determine 
an effective temperature in the range of 43000K to 48000K. Smith et al. favored the lower 
range for IRS2, 43000 — 44000 K, corresponding to a spectral type of 05 — 05.5 (Panagia 
1973), but not ruling out a temperature of 48000 K, corresponding to an 04.5 spectral 
type. More recent modelling suggests the lines used by Smith are very sensitive to model 
parameters, and after an effective temperature of 40000K, the lines are rather insensitive 
(Mokiem et al. 2004). Therefore, a new look at Smith et al. suggests only that the effective 
temperature is greater than 40000K. 

The ionizing flux can be measured independently at radio wavelengths. At 5 GHz, 
Wilson et al. (1970) measured a flux of 172 Jy (over a large area with a 4' beam), while Vigil 
(2004) measured 99 Jy (with a 10" beam interferometrically) at the same frequency. The 
Vigil (2004) observations focused on the immediate region around IRS2, and is likely to be 
more representative of the flux solely due to IRS2. 

Taking the electron temperature to be 7500 K (Caswell & Haynes 1987), and assuming 
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the Wilson et al. and Vigil results represent the upper and lower limits respectively, we find 
that the number of Lyman continuum photons, Ni, is in the range 2.7 — 4.7 x 10^^ s~^ 
(Armand et al. 1996), or log {Nl) ~ 49.1 - 49.4. Following Martins et al. (2005) this range 
of log {Nl) values implies a spectral type of 04 to ~05.5 for luminosity class (LC) V stars, 
and 05.5 to 07 for LC III stars. For an effective temperature > 40000 K, the calibration of 
Martins et al. (2005) suggests a spectral type of 05.5 or earlier for LC V and 05 or earlier for 
LC III. Therefore, the radio fiux and effective temperature measurements give inconsistent 
results for LC III, while for LC V they are consistent. Futher, as these are young stars, LC 
V is more likely, and thus we conclude that the spectral types of the IRS2 binary, which we 
assume to be of equal mass, are likely to be ~05 V, with a range of 04-05.5 consistent with 
the available data. Note that the other OB star candidates in the immediate vicinity of IRS 
2 have log {L^oi) < 3.4 (Wolk et al. 2006), and thus are not likely to contribute significantly 
to the continuum radiation. 



4.3. Cluster Membership 

Completeness histograms of H and Kg magnitudes do not show the characteristic steep 
rise towards the fainter magnitudes that is often present due to background sources if the 
cloud is being penetrated. This implies that the region is relatively free from extensive back- 
ground contamination to the depth of our survey (Carpenter et al. 1997), but the likelihood 
of a small fraction of the detections being from background sources is estimated from the 
SOFI control field. 

An extinction map for the NACO field was calculated following the method of Guter- 
muth et al. (2005) and is shown in Figure O We used H — Kg colors for groups of the 20 
closest sources detected in all bands, sampled at each point in a uniform grid with a spacing 
of 3'.'7, in order to analyze the extinction surface density. The reddening law of Rieke and 
Lebofsky (1985) was used. Sources used are all sources detected in H Sz Kg with a < 0.2. It 
is important to note that the measurements used to make this extinction map are dominated 
by cluster members, and so it is mostly a map of the extinction in front of and intrinsic to 
the cluster. 

The NACO extinction map yields a mean extinction of Ay = 14.8, which is a much 
higher extinction than the SOFI field due to the molecular cloud. The average extinction 
in the SOFI Kg-hand control field is Ay = 1.83, which was calculated as described above 
for the NACO field. Therefore, for comparison, the control field Kg-hand magnitudes were 
artifically extinguished by Ak = 1.45 magnitudes. This is an underestimate of the reddening, 
thus the amount of contamination should be considered an upper limit. The extinguished 
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sample yields a 90% completeness limit of 17.12 magnitudes for the control field compared 
to 16.97 for the RCW 38 field. 

Figure [HI shows Kg magnitude histograms (regardless of a detection or lack of at J and/or 
H) for the NACO RCW 38 data, the artificially extinguished control field data scaled by 
area, and the "corrected" histogram given by their subtraction. The total FOV of the 
control field is 4'.2 x 4'.2 and the area covered in the NACO observations is only 5.8% of this, 
which gives the area scaling factor. Integrating over the bins brighter than the NACO 90% 
completeness limit cutoff and using ^/n counting statistics to calculate error, the histogram 
shows 424 total NACO K^-band sources and N^kgd = 20 ± 2 control field sources (scaled by 
the area ratio) so that 4.7% of the NACO sources are potential background contaminants. 
Therefore, of the 337 JHKg sources (including protostellar candidates), 317 ±5 sources can 
be considered cluster members. This estimate of background contamination should represent 
an overestimate because the number of reddened background stars is reduced by requiring 
J and H detections as well. Because the contamination from background sources is small, 
all 337 sources are included in the analysis and diagrams of this paper, except for the disk 
fraction analysis where the corrected total of 317 ± 5 sources is used. 

The small amount of contamination is due to two reasons. One is simply the small 
field-of-view of our observations, combined with the extinction in front of and intrinsic to 
the cluster, of Ay ~ 15. The other is the large background column density, as derived from 
observations of the thermal dust continuum near 1-mm wavelength (Cheung et al. 1980). 
The peak column density is measured to be 8 x 10^^ cm~^, corresponding to an extinction 
of Ay = 400 (Wolk et al. 2006). Thus it appears that the embedded cluster lies in front of 
a dense molecular cloud whose high extinction screens out background sources. 

4.4. Candidate Protostars and Kg-hand only Sources 

There are five JHKg sources whose la error bars in both coordinates of color-color space 
fall completely outside the reddening vector extending from the T-Tauri locus in Figure [21 
These are protostellar candidates identified in a similar way as Lada & Adams (1992), and 
they represent approximately 2% of the population. They are tagged in Table [3] with a 
'p'. These candidates were visually inspected and appear to be true detections. Lada et 
al. (2000) find that protostellar candidates compose 13% of sources detected in JHKL in 
the Trapezium Cluster, with the difference likely due to the additional L-band sensitivity 
represented in the Lada et al. data set. 

There are also 53 sources with magnitude error less than 0.2 that were detected in Kg- 
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band only. These sources are mostly cluster members due to the small amount of background 
contamination as shown in Figure [61 These sources could either be highly extinguished T- 
Tauri stars or protostellar candidates. We have compared their distribution to the extinction 
map shown in Fig \5[ and we find that most of the sources are located in regions of high 
extinction. In Figure [1] these sources are clearly seen in the east and south-east as red 
sources. There are also 65 sources with magnitude error less than 0.2 that were detected 
in H and Kg-hand only. These sources are distributed in a similar way to the sources only 
detected in the fCg-band, suggesting they are likely extinguished T-Tauri stars or embedded 
sources. 



4.5. Molecular Hydrogen Emission 

Narrowband imaging of the H2 1-0 S(l) line at 2.122 fim is often used to identify 
shocked molecular emission due to outfiows. In a nebulous region such as RCW 38, pure H2 
line emission must be distinguished from the extensive and bright continuum emission using 
multi-band data. By comparing the 2.122 fim narrow-band image with the broad-band Kg 
image, two regions of H2 emission have been identified (Fig. [7]). Both lie along the IRS 1 
ridge, which has previously been identified as a region of active star formation (Smith et al. 
1999; Wolk et al. 2006). Region A in Fig. [7] shows H2 emission extending to the SW from 
an infrared point source that is one of the redder sources in the region (No. 174 in Wolk 
et al. 2006). This infrared source is coincident with one of the brighter X-ray sources (over 
500 counts) but is barely detected at J in the NACO data. It is likely to be an embedded 
protostar with an outfiow. 

In region B, H2 emission extends to the SW from a kink in the IRS 1 ridge emission, 
suggestive of a relationship between the two features. No infrared point source can be 
directly associated with this emission however, and so its origin is unknown. A deeply 
embedded protostar may be waiting to be found in this region, perhaps through deep mid- 
infrared imaging and/or through interferometric imaging at millimeter wavelengths. Unlike 
the ONC (Kristensen et al. 2003; McCaughrean & Mac Low 1997; Allen & Burton 1993), or 
a large number of small scale outfiows from the cluster members (Davis et al. 2008), there 
is no large scale H2 emission in the RCW 38 region that could be suggestive of a recent 
explosive event, or a large number of small scale outfiows from the cluster members. 



4.6. Photoevaporating Disks 



The narrowband imaging reveals tliree candidate photoevaporating disks - stars with 
comet ary-like head-tail morphologies (Figs. [H [9], [TOl [TTj) . While all have stellar-like point 
sources at their heads, only one of these is clearly directly affected by IRS 2, with its tail 
pointing directly away from the 0-star and lying in close proximity (Fig. These objects 
are well seen in the narrowband images at 2.12 /zm (H2) and 2.17 fim (Br7), when compared 
to the broadband images at H and Kg that are dominated by continuum emission. The 
proplyds in Orion have typical sizes C.'IS to 1" (O'dell 1998), which translates to 62-413 
AU assuming the distance to the ONC is 414±7 pc (Menten et al. 2007). Our observations 
have a pixel scale of 0'.'054 (91 AU at a distance of 1.7 kpc), which should allow for the 
identification of any larger Orion-like proplyds in the region. We have detected at most 
three similar objects with tail sizes that are approximately 1275 AU (Fig. [TOj) and 2208 AU 
(Figs, m [11]). It is interesting to note that two of these objects are located in the IRS 1 ridge. 
One of the goals of undertaking both narrowband and broadband imaging was to identify 
both resolved and unresolved proplyds. In the case the proplyds are unresolved, they would 
be identified as sources of line emission in the narrowband images after careful subtraction 
of the continuum (broadband) emission. However, the variable PSF, both in time and across 
the FOV, make such a comparison difficult if not impossible, and we have not attempted to 
do this. Stable PSFs hke those provided by space-based platforms are needed to undertake 
a search for unresolved emission from evaporating disks. 



4.7. Disk Fraction 



Using the (J — H) vs. {H — Kg) color-color diagram (Figure [2]), it is possible to isolate 
87 sources (92 including protostellar candidates) with Kg-hand emission greater than that 
expected for a stellar photosphere. These sources are tagged in Table [3] with an 'e' for excess 
sources, and a 'p' for protostellar candidate. Figure [5] shows a spatial distribution of these 
sources and Figure [12] shows a histogram of sources vs. radial distance from IRS2. Excess 
infrared emission is often indicative of the presence of inner disks. The disk fraction D is 
calculated as: 



^excess 

Ntot — Nfjkgd 

where N^xcess is the number of sources with an infrared excess, Ntot is the total number 
of sources, and Nbkgd is the number of background sources. The disk fraction as a function of 
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radial distance from IRS2 is plotted in Figure [131 In this manner we find a disk fraction of 
29 ± 3%. The disk fraction is quoted without mass and extinction limited sampling because 
it did not affect the final value. We perform mass and extinction limited sampling when 
comparing RCW 38 and the ONC. This is a lower fraction than found in similar studies 
(JHK) of similar clusters (Lada et al. 2000; Haisch et al. 2000; Hillenbrand 2005), given its 
degree of embeddedness {Ay ~ 10) and inferred young age (0.5 Myr; Wolk et al. 2006). 

As mentioned previously, RCW 38 is an interesting region to study in that we are able 
to compare what we find to the best studied young cluster, the ONC, to see if the trends 
found there exist elsewhere. Lada et al. (2000) in their JHKL survey of the ONC covering 
approximately the same linear scale as our RCW 38 study (a 0.4 pc square region) found a 
disk fraction of 50 ± 7% determined by X^-band excess, while a disk fraction of 80 ± 7% was 
determined using L-band excess. Both Kg and L are sensitive to warm inner disk emission 
(< 1 AU) and require favorable angles of inclination, inner disk holes, and accretion rates 
for detection (Lada et al. 2000). However, the photospheric emission of a given star has 
a lower intensity in L rather than Kg., which means the X-band excess must compete with 
photospheric emission much more than L-band. Thus, longer wavelengths such as the L- 
band, that can more reliably detect disks, resulting in higher disk fractions than for studies 
relying on iTg-band excess only, are needed. It is highly likely that our survey of RCW 38 
only detects a fraction of disks present. The only way to test this would be to image the 
region in L-band or longer wavelengths with sufficient sensitivity and resolution to resolve 
individual sources, which is a difficult task for a cluster at 4 times the distance of the ONC 
(1.7 kpc) with bright nebulosity. 

The disk fraction of 50 ± 7% found for the ONC using Ks-haxid excess is difficult to 
directly compare to the present result due to differences in sensitivities and photometric 
techniques. In order to make a fairer comparison, the JHKg data on the ONC from Muench 
et al. (2002) was obtained and a 0.3 Mq mass cut and a (J — H) < 2 extinction cut applied 
to both data sets, over the same linear region. This resulted in 218 sources in the ONC, 
including 71 excess sources, and 112 sources in RCW 38, with 28 excess sources. The overall 
disk fractions for the two data sets subjected to these cuts are 25 ± 5% for RCW 38 and 
42 ± 8% for the ONC. These fractions are significantly different, for reasons unknown. As 
noted, however, determining the infrared excess fraction from JHK data only is likely to 
underestimate the true fraction (Lada et al. 2000). 
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5. Analysis 

5.1. Photoevaporating Disks? 

Sources with disks in RCW 38 seem to be evenly distributed throughout the region, with 
the majority of disks found between 0.07 and 0.17 parsecs from 1RS2 (48%, compared with 
only 24% at smaller distances). Because of the UV radiation emanating from the massive 
0-star binary IRS2, it might be expected that the disks nearest to the central source would 
have been evaporated. Photoevaporating disks in the inner region of a cluster have been well 
documented in the Trapezium Cluster on account of the O star 9^ C Ori through both radio 
observations (Churchwell et al. 1987; Felh et al. 1993) and both ground and space based 
optical imaging (Laques & Vidal 1979; O'dell, Wen, & Hu 1993). Instead of a deficit of disks 
near IRS2, we find a slight increase in disk fraction with decreasing projected distance from 
IRS2 (Fig [T3l) . However, the physical separations are not known and the proximity of some 
of these stars to IRS2 could be a projection effect. In addition, since the cluster is quite 
young and thus dynamically active, stars moving at 1 km/s would move across the entire 
FOV in just 0.5 Myr, making it nearly impossible to tell which stars have spent prolonged 
amounts of time near IRS2. This could explain the apparent lack of resolved proplyd-like 
objects in the inner region near to IRS2. 

It is also possible that the outer disks of many of these sources have been blown away, 
leaving the inner disks intact which is the portion most sensitive to Kg-band measure- 
ments. 05.5 stars such as IRS2 emit large quantities of Far-Ultraviolet (FUV) and Extreme- 
Ultraviolet (EUV). More FUV radiation is emitted than EUV radiation, which means that 
FUV radiation is the dominant cause of photoevaporation because it works to evaporate the 
outer portions of the disks at a faster rate. EUV radiation works to evaporate the disk at a 
slower rate, so when the outer portion of the disk has been photoevaporated from FUV radi- 
ation, it is the more energetic EUV radiation that slowly works to cook away the inner disk. 
Adams et al. (2006) modeled this photoevaporation effect for embedded clusters with varying 
memberships. They define a Go term given in units of the typical interstellar radiation field 
at FUV wavelengths (given as 1.6 x 10~^ ergs s~^ cm~^). For an O star binary (assuming a 
temperature of 41000K and i?* = 15i?0), sources at 0.01 pc would experience Gq = 43000 
which decreases as l/(distance)^ so that sources at 0.1 pc would experience Gq = 430. Since 
the O stars represent the dominant source of FUV radiation, the Gq values given take into 
account just the central source and ignore contributions from nearby less massive sources. 
However, if we were to take into account the 31 candidate OB stars identified in Wolk et 
al. (2006), the average FUV radiation experienced by a given star would be greater. Higher 
Go values represent a faster photoevaporating effect and most easily affect lower mass stars 
whose disks are less tightly bound. 
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If we assume large planet formation occurs in the region of 5-30 AU, for M^, = 0.5Mq 
or less it would be inhibited over a 2 Myr period in a region where Gq = 3000, though this 
would only impact the stars minimally if the cluster is < 2 Myr. As an example from Adams 
et al. (2004), disks around 0.25 Mq stars would be photoevaporated down to around 14 AU, 
reducing the reservoir of material for planet formation around the star by half and likely 
preventing the formation of outer planets. Without velocity dispersion data for the cluster, 
however, it is difficult to know what kind of orbits the cluster members follow around IRS2 
and how much time they actually spend close to this central source. However, if we assume 
circular orbits for cluster members so that they spend 100% of their time at their projected 
radial distance, there are 18 JHKg sources that experience Gq > 3000. Of these, 5 sources 
have detected infrared excess, giving a disk fraction of 27%. However, random sampling of 
the sources using the disk fraction calculated previously gives a 28% chance of choosing 5 
or less excess sources, so that this slightly lower disk fraction is not statistically significant. 
This result is to be expected since i^'s-band excess detects the inner wall of disks which can 
be found at distances less than 1 AU from the central source so that even if a star's disk 
was significantly reduced in size, this lack of disk material would not be noticeable. As an 
example, Eisner et al. (2008) measured the mass of disks in the ONC versus Taurus and 
found that the ONC disks are much less massive, indicating that photoevaporation effects 
from the OB stars in the field, have played a key role in evaporating the outer disks in the 
ONC. 



5.2. Cluster Structure 

5.2.1. Azimuthally Averaged Surface Density Profile 

The azimuthally averaged surface density profile, or radial profile, of a cluster is a 
common tool used to probe the decrease in surface density as a function of radius from the 
central peak (Carpenter et al. 1997; Hillenbrand & Hartmann 1998). King (1962) surveyed 
15 globular clusters in order to examine their radial profiles and found a relationship that 
could be expressed with only the tidal radius and the core radius of a cluster. Hillenbrand & 
Hartmann (1998) fitted their radial profile for the ONC using such a King profile. Neither 
the ONC or RCW 38 is a globular cluster in equilibrium which makes the King model less 
physically significant for this application; however, as RCW 38 shares many similarities with 
the ONC, it seemed reasonable to determine the King profile parameters for RCW 38 for 
comparison to the Hillenbrand & Hartmann values. 

The radial profile of RCW 38 was made using only the sources detected at Ks with 
magnitude < 17 and uncertainties < 0.2. As shown in Figure [H] it is relatively fiat inside 
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the core of the cluster at radial distances less than 0.1 pc and then begins to fall off at further 
distances. The radial profile can only be examined up to a distance of 0.25 pc from the center 
of the cluster in order to limit edge effects caused by sampling distances outside the FOV 
in certain radial directions. This radial profile represents only the sources resolved in this 
sample, whereas it is possible that there are deeply embedded and extinguished sources that 
are cluster members but not taken into account in this radial profile. Figure [15] shows a 
contour map of Kg stellar surface density that illustrates the fairly uniform distribution with 
only a few small areas of possible sub clustering. 

The shape of the radial profile in Figure [H] motivates the choice of a flat-topped profile 
with a power-law drop off to fit the radial profile of the data such as the King model. For 
regions where the tidal radius is much larger than the core radius (as is generally assumed), 
the King model can be expressed as: 



where tq is the core radius of the cluster, and ctq is the initial stellar surface density. 
With Tc = 0.1 pc and ctq = 5700 pc~^, the model provides a reasonable fit to the radial profile 
of RCW 38 (Fig. [T4|) . Hillenbrand & Hartmann (1998) give best-fit King model parameters 
of To = 0.164 pc and cxo = 5600 pc~^ for the ONC. This central density is very close to 
that of RCW 38 while the scale radius is larger, suggesting that RCW 38 is more centrally 
concentrated than the ONC. 

Because the King model is applicable for globular clusters, observations show that it is 
not a good fit to the extended outer envelopes of young star clusters in our Galaxy and the 
LMC and SMC. In this case a modified King model (Elson, Fall, & Freeman 1987; EFF) is 
a better description: 



where the introduction of the exponent 7 allows the more extended nature of Galactic 
clusters to be fitted. This model differs from the King model in that there is no truncation 
radius included. For RCW 38 we obtain a best fit with = 0.1 pc, do = 5200 pc~^ and 
7 = 1.7. As in the case of the King model, the EFF model provides a reasonable fit to the 
data in the fiat and power-law fall-off regions. 




1 



(3) 




(4) 
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5.2.2. Nearest Neighbor Distance Analysis 

We use a nearest neighbor distance analysis to look for subclustering among close pairs. 
Our King model (eqn. [3l dashed line in Figure [H]) was used to generate a group of 1000 
random clusters with the same radial distribution and number of sources as ROW 38. If 
the RCW 38 distribution of close pairs was significantly different from the Poisson noise 
models, then we would expect subclustering among close pairs. A nearest neighbor distance 
histogram of the NACO data with a bin size of 800 AU is compared with the averaged 
histogram values for these model clusters in Figure [161 The model fits with the data quite 
well, with almost every bin matching when taking into account the error bars. This shows 
that there is a fairly low degree of subclustering on a smaller scale, a fact that should be 
expected given the fairly flat radial profile. 

5.2.3. Q Parameter Analysis 

Cartwright & Whitworth (2004) use a Q parameter in order to classify star clusters 
based on the degree a cluster is dominated by either fractal subclustering or a monolithic 
radial density gradient. They show that a uniform 2D (or unconstrained 3D projected into 
2D) distribution gives Q ~0.72 while a uniform 3D distribution constrained to a spherical 
volume and projected into 2D gives Q ~0.79. Thus, spherical clusters with some degree of 
subclustering are expected to have Q < 0.79 while clusters with a smooth radial gradient 
are expected to have Q > 0.79. 

The parameter Q is defined as Q = ^ where m is the mean distance from a minimum 
spanning tree scaled by a factor of (with N being the number of stars and A being the 
2D cluster area) and s is the mean of the two point correlation function measured for all of 
the stars in the cluster divided by the cluster radius. Using this procedure in order to further 
characterize the structure of RCW 38, a value of Q =0.84 is found, which places it on the 
border between clusters that are characterized by a relatively uniform distribution (smooth 
with a large scale radial gradient) and clusters that show some degree of subclustering. 
Together with the analysis in the previous section, this result suggests that RCW 38 is quite 
uniform and does not show strong subclustering. 

5.2.4- Larson Analysis 

Figure [T71 recreates the surface density vs. angular separation plot from Larson (1995). 
The power law for RCW 38 differs fairly significantly from the original Larson plot and does 
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not contain the same elbow evident in the Larson plot that is identified as the transition 
from the clustering to the binary regime. The absence of the elbow is easily explained by the 
fact that it occurs at an angular separation of 8663 AU, which is at the very high end of the 
nearest-neighbor histogram in Figure [in] and thus not easily measurable by this smaller field 
of view. Bate et al. (1998) also explains that we should expect the elbow in the Larson plot 
to be washed out in a large 3D cluster where the angular separation of many sources can 
be attributed to projection effects. In addition, from Figure [17] it is clear that RCW 38 is 
much denser than the Taurus cluster and thus it is possible that different physical processes 
govern star formation in these regions. 

Nakajima et al. (1998) uses a similar Larson analysis on five young star forming regions 
observed at visible wavelengths, fitting broken power laws to their sample. As a similar break 
in the Larson plot is not observed in our RCW 38 data, due to the high density, we can 
only compare the slope found at larger angular separations. Table H] shows the slopes from 
Nakajima et al. (1998), where 75- gives the power-law index at small angular separations 
while 7l gives the power-law index at large angular separations, alongside the Larson result 
for Taurus and the value for RCW 38. The RCW 38 result is intermediate between the 
values of 75- and 7/, listed for other regions. 

It is difficult to compare the deep near-infrared data from RCW 38 to the V and R band 
data presented by Nakajima et al. because in those bands there is poor completeness due 
to the high levels of extinction in the visible band [Ay ~ 10-30). Because of the variability 
in extinction and distance among the clusters surveyed, the mass sensitivities in each region 
are likely to differ significantly. Allen et al. (2002) show that these factors can affect the 
Larson plot results significantly in their HST/NICMOS survey of the Ophiuchus cluster. 

6. Summary 

We have presented high resolution adaptive optics near-infrared imaging (broadband J, 
H, Kg] narrowband 1.28 /xm, 2.12 fim, and 2.17 ^m) of the deeply embedded massive cluster 
RCW 38, over a field-of-view of ~ 0.5 pc, centered on the 05.5 star IRS2. Our main results 
are: 

1. A total of 360 sources are detected in all three broadbands. Of these, 344 have uncer- 
tainties in all bands of < 0.2 mag, and the remainder are removed from the sample. A 
further seven sources were removed due to contamination from other brighter nearby 
stars or due to image artifacts, leaving 337 reliable sources. Of these 337, five have 
H — Kg colors in excess of that expected for T Tauri stars, and are considered to be 
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protostellar candidates. A comparison to a nearby control field suggests that 15-25 
sources are unassociated with the chistcr. There are 53 Kg only detected sources, that 
are most likely to be cluster members due to the small amount of background contami- 
nation. These sources could either be highly extinguished T-Tauri stars or protostellar 
candidates. 

2. A J — if vs. H — Kg color-color diagram is used to identified stars with an infrared 
excess indicative of inner disk emission. In this manner the disk fraction is found to 
be 29 ± 3%. For comparison to similar data from the ONC, a 0.3 M0 mass cut and 
a (J — H) < 2 extinction cut was applied to both regions over the same linear area, 
resulting in disk fractions of 25 ± 5% for RCW 38 and 42 ± 8% for the ONC. 

3. The central star IRS2 is found to be an equal mass binary of spectral type ~05 and 
projected separation 459 AU. 

4. Molecular hydrogen line emission, a good tracer of shocked outflow emission, is seen 
at two locations, from a comparison of the narrowband 2.12 /xm images with both 
broadband Kg and narrowband 2.17 /im images. At one location this emission has a 
morphology suggestive of outflow emission, and is associated with a deeply embedded 
source that is also a hard X-ray source (Wolk et al. 2006). 

5. Three candidate photoevaporating disks are identified, with head-tail morphologies, 
and with stellar-like sources at their heads. Only one lies in close proximity to IRS2 
and has a tail pointing directly away from the O star. Although one goal of the study 
was to identify unresolved proplyds through a careful comparison of narrowband and 
broadband imaging, this was not possible due to the varying PSF between the bands. 

6. The cluster radial density profile is well described by a modified King profile, which is 
widely used to fit globular clusters. In young star clusters, such a profile breaks down 
in fitting the extended outer envelopes. This is not a problem for the data presented 
here, as the FOV is too small to trace the outer envelope region (~ 0.25 pc). Nearest 
neighbor analysis, both directly on the data and on a model density profile that fits the 
data, show essentially no substructure (sub-clustering). Comparison to other clusters 
suggests that this is not surprising for our limit FOV. The stellar surface density in 
the inner region of the clusters is ~ 5700 pc~^. 

There are a number of observational programs on RCW 38 that need to be undertaken 
in order to improve on this work and that of Wolk et al. (2006). The comparatively low 
disk fraction of ~ 29% found using only JHKg data needs to be confirmed through sensitive 
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observations in the L-band. The large distance and nebula emission will make these obser- 
vations difficult. The nature of the protostellar candidates and Kg-hand only sources should 
be further investigated, through sensitive mid-infrared imaging. The OB star candidates 
identified in Wolk et al. (2006) should be observed with infrared spectroscopy to determine 
their spectral types. Finally, the age of the cluster needs to be better determined, a difficult 
task given its extinction and youth. 

We thank Joana Ascenso for providing the SOFl data in advance of publication, and for 
sharing her ISAAC results in advance of publication. We thank the VLT Science Operations 
team for successfully executing our NACO observations in service mode. Partial support 
for this work was provided by NASA through contract 1279160 issued by JPL/Caltech, and 
by NASA contract NAS8-03060. This research has made use of NASA's Astrophysics Data 
System. 
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Table 1. DITs (Detector Integration Times) and NDITs for observations in each band. 



Band 


DIT (s) 


NDIT 


J 


4.0 


16 


H 


0.75 


80 


Ks 


0.50 


120 


NB_2.12 


10.0 


14 


NB_2.17 


10.0 


12 


NB_1.28 


90.0 


4 
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Table 2. Number of detected sources. 



Band 


Detections with 0.2 magnitude uncertainty cap 


J 


355 


H 


420 


Ks 


465 


JHKs All 


344 
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Table 3. Infrared Photometry of sources detected with magnitude uncertainties < 0.2. 



Source Number R.A. (J2000) Dec. (J2000) J J^rr H H^,, K, {Ks)err Notes 



1 


8 


59 


2.49 


-47 


30 


20.7 










14.09 


0.06 




2 


8 


59 


2.50 


-47 


30 


52.6 










15.02 


0.09 




3 


8 


59 


2.55 


-47 


30 


39.9 










14.41 


0.06 




4 


8 


59 


2.56 


-47 


30 


46.2 










12.08 


0.02 




5 


8 


59 


2.57 


-47 


30 


53.5 










13.09 


0.03 




6 


8 


59 


2.60 


-47 


30 


31.7 










13.35 


0.04 




7 


8 


59 


2.62 


-47 


30 


35.2 










12.86 


0.03 




8 


8 


59 


2.63 


-47 


30 


28.7 










15.07 


0.08 




9 


8 


59 


2.66 


-47 


30 


42.3 










13.52 


0.04 




10 


8 


59 


2.67 


-47 


30 


57.4 










15.68 


0.14 




11 


8 


59 


2.67 


-47 


30 


21.1 










14.83 


0.08 




12 


8 


59 


2 74 


-47 


30 


43 4 


17 72 


07 






1 5 32 


09 




13 


8 


59 


2.76 


-47 


30 


53.8 


17.28 


0.13 






14.95 


0.07 




14 


8 


59 


2.78 


-47 


30 


46.8 


14.97 


0.02 


13.61 


0.02 


12.65 


0.03 


e 


15 


8 


59 


2.83 


-47 


30 


34.0 










16.49 


0.15 




16 


8 


59 


2.89 


-47 


31 


6.5 










15.25 


0.10 




17 


8 


59 


2.91 


-47 


30 


30.8 


16.97 


0.05 


15.11 


0.05 


13.83 


0.04 


e 


18 


8 


59 


2.92 


-47 


31 


2.1 


15.01 


0.02 


13.44 


0.02 


12.36 


0.02 


e 


19 


8 


59 


2.93 


-47 


30 


54.1 


12.27 


0.01 


10.51 


0.01 


9.81 


0.01 




20 


8 


59 


2.93 


-47 


30 


54.1 


12.14 


0.02 












21 


8 


59 


2.97 


-47 


30 


39.8 


15.10 


0.02 


13.71 


0.02 


12.85 


0.03 


e 


22 


8 


59 


2.98 


-47 


30 


51.7 


17.88 


0.10 


15.90 


0.06 


14.12 


0.05 


P 


23 


8 


59 


2.99 


-47 


31 


0.9 






16.78 


0.10 


15.49 


0.10 




24 


8 


59 


3.00 


-47 


30 


44.4 


16.25 


0.03 


13.54 


0.02 


11.00 


0.01 


P 


25 


8 


59 


3.02 


-47 


30 


35.6 


14.88 


0.02 


12.67 


0.01 


11.47 


0.01 





Note. — In the Notes column, 'c' refers to those sources with an infrared excess, 'p' to protostars, 
which are sources with an excess greater than of that expected of a reddened T Tauri stars, and 'x' 
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to sources that have J — H < 1 and whose photometry is incorrect, for reasons described in the 
text. Units of right ascension are hours, minutes and seconds, and unites of dechnation are degrees, 
arcminutes, and arcseconds. Units of J, Jerr H, i/err Ks (i^s)err are magnitudes. Table[3]is pubhshed in 
its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance 
regarding its form and content. 
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Table 4. Slopes for Larson plots for Taurus (Larson 1995) and the clusters surveyed by 
Nakajima et al. (1998). 75 gives the power-law index at shorter angular separations while 
7l gives the power-law index at longer angular separations. 



Region 


7s 




Orion OB 


-1.6±0.4 


-0.15±0.02 


Orion A 




-0.23±0.02 


Orion B 




-0.69±0.01 


p Oph 


-2.5±0.3 


-0.36±0.06 


Chal 


-2.1±0.2 


-0.57±0.04 


Cha 


-2.4±0.5 


-0.55±0.03 


Vela 




-0.61±0.02 


Lupus 


-2.2±0.4 


-0.82±0.13 


Taurus 


-2.15 


-0.62 


ROW 38 




-1.32 
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Fig. 1. — Color composite image of RCW 38 through J (blue), H (green), and Kg (red) 
filters. The bright star in the center is IRS2. The blue sources to the north and south of 
IRS2 are J-band artifacts, as is the ring-like structure to the west of IRS2 in the nebulous 
filament (known as IRSl from mid-infrared wavelengths). 
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Fig. 2. — Color-Color Diagram for the NACO field. Plotted are the 337 sources detected 
in JHKg as discussed in Section 14.11 Reddened sources with a < 0.2 are displayed with 
crosses. Open circles indicate Kg-haiid excess sources identified by discarding sources within 
Icr of the lowest reddening vector in either dimension. Filled triangles represent protostellar 
candidates (Class I sources) or extreme CTTS's (Lada & Adams 1992 or Carpenter et al. 
1997). Upper and lower reddening vectors are from Rieke & Lebofsky (1985) and plotted as 
solid lines for an M5 giant (Bessell & Brett 1988) and M5 dwarf (B. Patten 2004, private 
communication) respectively. Reddening of Ak = 1 and 2 is indicated by the filled circles. 
The classical T-Tauri stars (CTTS) locus from Meyer et al. (1997) is plotted as a dotted 
line. Some of the sources are identified as non-excess sources despite appearing in quite far 
into the region where they would appear to be excess sources, due to greater than average 
errors that overlap with the non-excess region to the left of the first reddening vector. These 
sources are all faint sources, and are evenly distributed at different radii from IRS2. 
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Fig. 3. — Color-Magnitude diagram for NACO field. Symbols as in Figure EJ reddened 
JHKs sources are crosses, Kg-band excess sources are circles, and candidate protostars are 
triangles. The single black diamond indicates the only source detected in J and H that 
was not detected in Kg (thus our sample is not biased towards detected excess in Ks). 
Completeness limits for J and H are represented as dashed lines. The curved line at left is 
the unreddened 1 Myr pre- main- sequence isochrone from Baraffe et al. (2002) with reddening 
vectors (Rieke & Lebofsky 1985) at 0.025 Mq, 0.08 Mq, and 1.0 M0 shown as downward 
sloping lines. From this diagram it is clear that the cluster is reddened by around Ak = 1- 
Despite the surrounding differential reddening that exists due to the nebulosity, the cluster 
is surrounded by a bubble of nebulosity with the central region evacuated due to the wind 
from the O stars. This causes the overall reddening to be fairly constant over the sources 
within the cluster core. 
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Fig. 4. — a) Magnified image of fRS2 from tlie narrowband 2.12 fim image. Contours are 
plotted at 60, 70, 80, 90, 95 and 100% of tlie peak flux in arbitrary units, b) Magnified 
image of IRS2 from the narrowband 2.17 fj,m image. Contour levels are the same as in a). 
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Fig. 5. — Shaded extinction contour map for the NACO FOV. White areas indicate areas of 
higher extinction, expressed in a hnear scahng where black represents Ay = 10 and white 
represents Ay = 20. Contours are given for A Ay = 1, starting at Ay = 10. i^<j-band excess 
sources and potential protostellar candidates as identified in Figure [2] are overplotted using 
the same symbols as that figure (white circles and white triangles respectively). The bright 
central source IRS2 is overplotted as a white star. 
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Fig. 6. — a) Histogram of NACO Kg-hand magnitudes divided into bins of 1 magnitude. 
All sources detected at Kg, regardless of detectability at J or H, are included. The dotted 
line represents the 90% completeness limit, b) Histogram of SOFT control field Ks-hand 
magnitudes divided in bins of 1 magnitude and normalized by the ratio between the NACO 
FOV and the SOFI FOV. The dotted fine represents the 90% completeness limit. The 
histogram has been shifted to the right by the amount SAk as described in the text, c) 
NACO histogram corrected by subtracting the SOFI control field histogram. Only bins 
within the 90% differential completeness limit for the NACO measurements are included. 




Fig. 7. — Molecular hydrogen emission in the 1-0 S(l) line is revealed in these images, which 
combine narrow-band images at 2.12 /im (red), 2.17 /im (green), and 1.28 /im (blue). In the 
left panel is shown the IRS 1 ridge region, and on the right are shown zoomed in regions A 
and B. The H2 emission is shown in red. 
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Fig. 8. — Finding chart for proplyd-like objects. The image is the narrowband 2.17 /xm 
image, and the three proplyd-hke objects shown in more detail in Figs. [9lfTT] are indicated 
by the boxes marked C, D and E. IRS2 is indicated with an asterix. 
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Fig. 9. — Image at 2.12 /iin showing a candidate proplyd with a tail near IRS2. Identified as 
Region C in Figure [H The arrow points to the head of this globule, with the tail stretching 
leftwards and slight upwards, away from the location of IRS2. 
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Fig. 10. — Images in four bands (as noted) of a proplyd candidate with tail. Identified as 
Region D in Figure [H 
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Fig. 11. — Images in four bands (as noted) of a proplyd candidate with tail. Identified as 
Region E in Figure [H The arrows point to the head of the globule. 
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Fig. 12. — Radial surface density plot IRS2. Bin size is 4". The dashed line indicates all 
JHKs sources and the solid line is iTg-band excess sources from Figure [21 
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Fig. 13. — Disk fraction vs. radial distance from IRS2. Bin size is about 4". 
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Fig. 14. — Radial stellar surface density plot for RCW 38 centered on IRS2 using NACO 
data. Bin size is about 0.02 pc (2'.'5). The stellar surface density is computed for all sources 
detected with K < 17 with a < 0.2. The dashed line represents the fitted King model 
(eqn. [3]) while the dotted line represents the fitted Elson et al. (1987; EFF) model (eqn. H]). 




Fig. 15. — Stellar surface density plot. The contours of Kg stellar surface densities for 
sources with K < 17 and a < 0.2 are shown in the upper image. The contours represent 
stellar surface densities of 30, 100, 500, 1000, 2500, 5000, 10000, 20000, and 25000 stars 
pc~^. Densities are computed using 5 nearest neighbors. For reference, in the lower panel is 
shown the same FOV using the Kg band image. 
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Fig. 16. — Nearest neighbor histogram. Black represents NACO data (using sources with 
K < n and a < 0.2) and grey is the average of 1000 randomly generated clusters with 
the same radial profile as shown in Figure [TH Error bars represent the standard deviation 
of the random averaged clusters in each bin. Bin size is 800 AU. As a note, the two-point 
correlation function (TPCF) of this same data set yields nearly identical results, so it is not 
included here. Since the data does not deviate significantly from the model (where sources 
are randomly distributed for the given radial profile), it is consistent with the model that 
no major substructuring exists in this FOV. This is as expected because the radial density 
plot in Figure [H] is essentially flat over most of the region. 
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Fig. 17. — Surface density of companions as a function of angular separation (after Larson 
1995). The Larson result, scaled to the same distance, is shown as a dashed line (the original 
Larson result shows an elbow at around log 9 = -1.75 which translates to log 9 = -2.83, for 
the distance to RCW 38 of 1.7 kpc, cf. 140 pc for Taurus). The grey line shows a power law 
of -1.32 which is significantly different from the -2.15 power law for the binary regime and 
-0.62 power law for the clustering regime measured in Larson (1995). 



